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The thermodynamic properties of non interacting bosons on a complex network can be strongly 
affected by topological inhomogeneities. The latter give rise to anomalies in the density of states 
that can induce Bose-Einstein condensation in low dimensional systems also in absence of external 
confining potentials. The anomalies consist in energy regions composed of an infinite number of 
states with vanishing weight in the thermodynamic limit. We present a rigorous result providing 
the general conditions for the occurrence of Bose-Einstein condensation on complex networks in 
presence of anomalous spectral regions in the density of states. We present results on spectral 
properties for a wide class of graphs where the theorem applies. We study in detail an explicit 
geometrical realization, the comb lattice, which embodies all the relevant features of this effect and 
which can be experimentally implemented as an array of Josephson Junctions. 

I. INTRODUCTION 

The impressive experimental demonstration of Bose-Einstein condensation (BEC) has stimulated a wealth of the- 
oretical work aimed at a better understanding of its basic mechanism and of its potential consequences for the 
engineering of quantum devices . 

A well known general result |Q| is that for an ideal gas of Bose particle, BEC does not occur for homogeneous 
systems in dimension d < 2; the same is true for non interacting bosons on regular periodic lattices. The result 
cannot be extended to more general discrete structures lacking translational invariance, which, in principle, can now 
be implemented by combining quantum devices in non conventional discrete geometrical settings. 

The first important problem one is faced with is the characterization of the geometrical properties of a general 
network, i. e. a graph. On one hand, one is interested in describing the complex geometry of a sample and identifying 
its effects on physics. On the other, once the effects of geometry are known, one can use the geometrical setting 
to tune the physical properties, by a sort of engineering in network design. As we will see, this can lead to effects 
analogous to the introduction of an external confining potential for regular geometries. 

The problem of describing the large scale geometry of graphs by an effective parameter that generalize the Euclidean 
dimension of regular lattices has been successfully solved by the introduction of the spectral dimension d j^Q] , which 
can be experimentally measured |^] and rigorously defined by graph theory ||. However, one expects the influence 
of topology to be richer and more complex on discrete structures, due to the possible relevance of local geometrical 
details, in addition to the large scale structure described by dimensionality. 

In this direction, recent works on BEC on inhomogeneous networks have put into evidence that strong inhomo- 
geneities can give rise to condensation at finite temperature for a system of non interacting bosons in absence of an 
external confining potential even in a low-dimensional structure as a comb lattice, where d = 1. This phenomenon 
arises from a peculiar part of the density of states on inhomogeneous networks we shall refer to as hidden spectrum 
||. This hidden spectrum consists of an energy region filled by a finite or infinite number of states which do not 
contribute to the normalized density in the thermodynamic limit. Hidden spectra do not usually affect bulk thermo- 
dynamic quantities but, as we shall show in the following, can have dramatic effects since bosonic statistics allows the 
macroscopic occupation of a single quantum state. 

In this paper we shall determine the conditions under which, even for low dimensional systems, non-interacting 
bosons on a general discrete network may lead to BEC ||. We shall see that this is indeed possible if one resorts to 
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a suitable discrete inhomogeneous ambient space on which bosons are defined; this is very desirable in view of the 
engineering of quantum devices. 

Bosonic models over generic discrete structures can be made experimentally accessible through the realization of 
arrays of Josephson junctions (JJA) |g|Jlo[|. The latter are devices that can be engineered to realize a variety of non- 
homogeneous patterns. We shall see indeed that classical JJA arranged in a non-homogeneous network provide 
a relevant example of the proposed mechanism for BEC, leading to a spatial condensation of bosons in a single state. 

The model Hamiltonian describing the low-temperature JJA is the Bose- Hubbard (BH) model On a generic 
graph, the Hamiltonian is given by 

H BH = Uj2 n2 l +J2 A « ( Vn * n 3 ~ J ( a l a J + a} *)) (!) 

i ij 

where Aij is the adjacency matrix of the graph: = 1 if the sites i and j are nearest neighbours and Aij = 
otherwise; a\ creates a boson at site i and = a\ai. The phase diagram reflects the competition between the boson 
kinetic energy (hopping term) and the repulsive Coulomb interaction. In a realistic experimental setup ]l2|| , the 
parameters U and V depend on the ratio between the intergrain capacitance C and the gate capacitance Co, while 
the parameter J describes Cooper pair hopping. Classical JJA are obtained when U, V « J (i.e. C/Co — > 0); in 
this limit the hopping term dominates the physics of the system. 

Recently, the BH model has been shown to describe the dynamics of ultracold bosonic atoms on a particular 
regular discrete geometry, the optical lattice p3| . In this context the BH model describes the hopping of bosonic 
atoms between low vibrational states of the optical lattice; the parameters U and J depend on the power of the laser 
beams used to realize the lattice. The dynamics of a BEC in such a kind of periodic potential can thus be described 
by a discrete nonlinear Schrodinger equation [H]. Also for atomic Bose systems, the experiments have reached an 
high level of accuracy with easily tunable external parameters and accurately tailored trapping profiles Jl5| ] . 

In this paper we will consider a model of pure hopping of non interacting bosons on general non regular geometries, 
such as complex networks. In particular, we will analyse the possibility of having a Bose-Einstein condensate at 
finite temperature, providing a general condition for BEC on a general network We will describe in detail some 
geometrical settings where inhomogeneities act as a confining potential, giving rise to BEC even in low dimensions. 

The extension of this analysis to the full BH model on complex networks is not trivial and it is still the object of 
on-going investigations. 

The paper is organized as follows: in section II we introduce the mathematical notations and some basic concepts in 
algebraic graph theory; in section III we define the pure hopping model and prove a general theorem giving the most 
general condition for the existence of BEC at finite temperature for non interacting bosons on a complex network. We 
discuss in detail the presence of "hidden" regions in the density of states in the thermodynamic limit; in section IV we 
apply the theorem to the case of the comb graph showing that BEC at finite temperature is present in 1-dimensional 
structures as a consequence of the effects of the inhomogeneous geometry on the density of states; finally in section 
V we consider a wide class of networks where the relation between BEC and complex geometry can be analyzed in 
detail. 



II. COMPLEX NETWORKS AND GRAPHS 



On complex networks it is not possible to introduce important mathematical tools typical of Euclidean lattices, 
such as the Fourier transform and the reciprocal lattice. Nevertheless general and effective mathematical techniques 
are provided by graph theory and the algebraic approach to graph topology Jl6| . In this section we briefly recall some 
basic definitions and notations. 

A graph G is a countable set V of vertices i connected pairwise by a set E of unoriented links (i,j) = In 
general sites represent degrees of freedom and links represent interactions between them. In particular, for bosonic 
models, vertices can be regarded as Josephson junctions or as sites in the optical networks, while edges represent the 
hopping probabilities between sites. 

A subgraph S of G is a graph whose set of vertices V C V and whose set of links E' C E. A path in G is a sequence 
of consecutive links {(?, fc)(fc, h) . . . (n, m)(m, j)} and a graph is said to be connected, if for any two points i, j G V 
there is always a path joining them. In the following we will consider connected graphs. Every connected graph G is 
endowed with an intrinsic metric generated by the chemical distance r< 3 -, which is defined as the number of links in 
the shortest path connecting vertices i and j. On a graph the generalized Van Hove sphere S 0-r C G of center o and 
radius r is the subgraph of G containing all i € G whose distance from o is < r and all the links of G joining them. 
We will call N ^ r the number of vertices contained in S _ r . The Van Hove spheres play an important role in the study 
of the thermodynamic limit on infinite graphs. The physical quantities are usually evaluated by first considering the 
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models restricted to the finite spheres and then taking the thermodynamic limit by letting r — > oo. The average in 
the thermodynamic limit <p of a bounded function defined on the vertices <fr(i) is: 



0S hm ^ . (2) 

r— too iV r 

With this definition, the measure \V'\ of a subset V of is the average value of its characteristic function Xv(*)> 
defined by XV CO = 1 if i € V and xv(£) = otherwise. It can be shown JlTj that for graphs with polynomial 
growth, i.e. when N Tt0 ~ r p for r — > oo, the thermodynamic limit is independent from the choice of the center of the 
sphere o. Since all realistic networks which can be embedded in a 3-dimensional space have polynomial growth, we 
restrict to this class graphs and we drop the index o of the center of the sphere. The graph topology is algebraically 
described by its adjacency matrix: 

A -I 1 if (3) 

On disordered structures, a useful generalization of the adjacency matrix can be considered by introducing a 
bounded distribution of hopping probabilities Uj which differs from link to link: 

V} 131 \ = if A i;j = [ 1 



with sup^js \tij \ < 00. 

The Laplacian matrix Ltj is defined by: 



— ZiSij A^ (5) 



where Zi = Y]j Aij, the number of nearest neighbours of i, is called the coordination number. We will consider graph 
with bounded coordination number, i.e. with max^ Zi < 00. 

Lij is the generalization to graphs of the usual Laplacian on a lattice, where Zj = z, Vz. L^j has some important 
spectral properties: its spectrum is real, non-negative and bounded. In particular is a simple eigenvalue of Lij 
corresponding to the constant eigenvector. Notice that, while on a regular lattice is diagonalized by the Fourier 
transform, this is not the case for a generic graph. 

A fundamental problem in graph theory is the characterization of the large scale geometry of the graph by the 
definition of a dimension which generalizes to infinite graphs the usual Euclidean dimension d of lattices. The 
spectrum of the Laplacian operator defines the spectral dimension d, which has been shown to be the extension of d 
in many physical problems. On a infinite graph, the spectral density of the Laplacian matrix is denoted by p(l). The 
spectral dimension J is defined by the asymptotic behaviour of p(l) at low eigenvalues ^Q: 

P(0 ~ i J/2_1 for ^0+ (6) 

Interestingly, d can be a real number, it can be experimentally measured ||, it depends only on the large scale 
topology of the graph Q and it is the direct extension to graphs of the Euclidean dimension in many physical 
problems. In particular, the long time asymptotic diffusion of classical particles Jp), the low frequencies spectrum of 
classical harmonic oscillators ||, the zero mass limit of the Gaussian model H,[| and the existence of spontaneous 
magnetization in continuous symmetry spin models |is}| are determined by the spectral dimension. 

III. BEC ON COMPLEX NETWORKS: THE GENERAL THEOREM 

The Hamiltonian for non interacting bosons on a graph G is: 

H = E h V a i a o ( 7 ) 
ijev 

where a\ and etj are the creation and annihilation operator at site i, with [aj,oj] = <%. The Hamiltonian matrix Ziy- 
is defined by: 
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hij = Uj + SijVi (8) 

The term ty, defined in (||), describes hopping between nearest neighbours sites. The diagonal term V$ takes into 
account a potential at site i and it satisfies a boundedness condition: sup \ Vi\ < oo. In particular, the theorem holds 
also when Vj = 0, i.e. when there is no on site potential. 

The thermodynamic limit on the graph is studied by first considering the finite Van Hove sphere S r and then letting 
r — > oo. This definition of the thermodynamic limit is very general and it applies to a general discrete structure. 
On graphs with symmetries (lattices, fractals, bundled structures, etc.) different shapes of Van Hove spheres can 
be introduced, preserving the symmetry of the network. It can be shown that for regular choices of the spheres the 
properties of the model for r — ► oo do not depend on the shape of the sphere. The Hamiltonian (Q) restricted to the 
sphere is: 

h s " = £ /<;», (9) 

i,jes r 

s s 
where h^ = hij if i and j belong to S r and h^ — otherwise. 

For each finite sphere of radius r let us consider the eigenvalues equation: 

J2 KMj) = E W) (io) 

jes r 

The normalized density of states p r (E) on S r is defined as: 

P r {E) = ^Y. 6 { E - E D ( n ) 

k 

where E r k are the N r eigenvalues of H^. 

We define p(E) to be the density of states of hij in the thermodynamic limit if: 

lim / \p r (E) - p(E)\dE = (12) 

Let E m = Inf (Supp(/?(£J))), where Supp(p(_E)) is the support of the distribution p(E). On a graph, the asymptotic 
behavior at low energies of the density of states is described by an exponent a: 

p(E) ~ (E - E^- 1 for E -» E m (13) 

The exponent a in general depends both on graph topology and on the external potential. Interestingly, when hij = 
and when the external potential Vi satisfy suitable regularity conditions, a coincides with the spectral dimension d. 
In particular this happens when Vi = 0. A hidden region of the spectrum is an energy interval [Ei,E 2 ] such that 
[Ei, E2] H Supp(jd(£J)) = and linv^oo N^ Ei E ^ > 0, where NT Ei E ^ is the number of eigenvalues of in the interval 
[Ei,E2\- Notice that in general NT Ei E ^ can diverge for r — * 00 and the eigenvalues can become dense in [Ei,E-i\ 
in the thermodynamic limit. Therefore this condition not only includes the trivial case of discrete spectrum but is 
far more general; an interesting example of this behaviour is found in the comb lattice without external potential Q| 
which will be studied in detail in the next section. 

We now define the lowest energy level for the sequence of densities p r (E), setting E 7 = lnik(Ek) and Eq = 
linv^oo Eq. In general, Eo < E m . If Eq < E m , then [Eq, E m ] is a hidden region of the spectrum, which will be called 
hidden low-energy spectrum. 

In the following we will consider models at fixed filling / = N/N r (N is the number of particles in the system). In 
the grand canonical ensemble the equation that determines the fugacity z in the thermodynamic limit is: 

f p r (E)dE 

f = i™ / ~ — 1 ~BE , ( 14 ) 



Setting Eq — we have that < z < 1. 

The integral in the equation ( |l4| ) can be divided into two sums, the first one considering the energies smaller than 
an arbitrary constant e and the second the energies larger than e: 

p r {E)dE 1 f p r(E)dE 

j-l e /3E_l 2^ z -l e (3El_ l + z -l e f)E _ X ■ 



fe=0 



E>e 
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We define: 



E k <c 



fc=0 

as the fraction of particles with energy smaller than e. Bose-Einstein condensation occurs in this systems if it 
exists a critical temperature T c > such that, for any T < T c , n e = lmv^oo n r c > k > 0, for all e > 0; i.e. 
uq = lim c _»o n e — k > 0. 

From definition (|l6|), no can be strictly positive only if linv^oo z(r) = 1. Indeed, if this limit is smaller than 1 
it follows that n e < (1 — z)^ 1 lim^oo N£ /N r , where N* is the number of state with energy smaller than e. If the 
ground state is not infinitely degenerate, lim e ^o hm r ^ 00 (7V e r '/ ; /V r ) = and then no = 0. 

Taking first the limit r — > oo and then e — ► in equation (Uq) we obtain: 

,. f (p r (E) - p(E))dE f p{E)dE 
f = n + hm hm / ' FX " + hm km / ' (17) 



Or- ►oo./ E>e Z~ 1 ef 3E — l t^Or^oo J E>t z~ 1 e l3E - 1 



Now, from the boundedness of (z 1 e^ E — 1) 1 for E > e and from the definition (p"2|), the first of the two limits in 
the right hand side of (|l7|) vanishes: 

f p(E)dE f p(E)dE 

where, again, no can be different from only if z = 1. 

The integral in equation ( |l8| ) is an increasing continuous function of z with < z < 1. If the limit: 

/ c (fl = km / -4^- (19) 



is equal to oo, tken z < 1 and no = 0. If tke limit is finite, / c (/3) is a decreasing function of (3 witk lim^oo / c (/3) = 
and lim^o fcifi) = oo. Tken, for a suitable /3 C , / C (/3 C ) = /. For (3 > /3 C , (i.e. T < T c ) z — 1 and n = / - j c (/3) > 
wkile for (3 < j3 c , (i.e. T > T c ) z < 1 and n = 0. 

From tke divergence or finiteness of tke limit (|l^) one obtains tke most general condition for tke occurrence of 
Bose-Einstein condensation on a grapk. 

First, if = Eq < E m (i.e. tke system presents a low-energy kidden spectrum) tke limit (|l^) is finite and tkere is 
Bose-Einstein condensation at finite temperature: 

W) < J P (E)dE = (20) 

On tke otker kand, wken E = E m tke value of tke limit ( |l9|) is determined by tke exponent a. Indeed if a > 2: 

f 5 CxEi^dE [ P (E)dE 

wkere 5 and c\ are suitable constants. Tkerefore in tkis case Bose-Einstein condensation occurs at finite temperature. 
For a < 2 we kave: 

r ,ns f S zcxE%- 1 dE 

f c ((3) > km / -i— = oo 22 

Jo pE+l — z 



and tkere is no Bose-Einstein condensation. Wken a — 2 we kave to consider tke logaritkmic correction to ( |13[ ) and 
it is possible to skow tkat tke limit (|l9|) diverges. 

Tke Hamiltonian (R) describes different models of non interacting bosons on grapks. Tke simplest example is tke 
discretization of tke usual Sckrodinger equation. In tkis case tke Hamiltonian is: /ly = #-Ly, wkere is tke 
Laplacian operator. Tkis is a pure topological model, completely defined by tke grapk geometry and it can be skown 
tkat Eq — E m , i.e. tkere are no kidden states in tke low energy region and a = d. Tkerefore tke occurrence of BEC 
is determined by tke spectral dimension (^), wkick kas been calculated for a wide class of discrete networks ]q, ^9|,pQ[ . 

On real condensed matter structures an interesting case is a pure bopping of non interacting bosons on grapks. 
Tkis kas been considered in M for tke description of tke Josepkson junction arrays in tke weak coupling limit witk 
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Hamiltonian matrix given by: hij — —tAij. In this case the behaviour of the model is much more complex and the 
presence of hidden spectra can give rise to Bose-Einstein condensation also in low dimensional structures @] . 

Finally, the most general application of our result concerns particles interacting with an external potential, described 
by the Hamiltonian: hij = —tAij + Vi^ij ■ Indeed, the introduction of a local potential strongly enhances the possibility 
of modifying the energy spectra in order to induce BEC in non regular geometries. 

In the next section, we will study in detail the pure hopping model defined on the comb graph, where the inhomo- 
geneous topology produces low energy hidden states which give rise to Bose-Einstein condensation in a 1-dimcnsional 
structure (a = d = 1). 



IV. PURE HOPPING MODELS AND BEC ON THE COMB GRAPH 



The comb lattice is an infinite graph, which can be obtained connecting to each site of a linear chain, called backbone, 
a 1-dimensional chain called finger (see figure 1). The sites of the comb can be naturally labelled introducing two 
integer indices (cc, y) with x, y £E where x labels the different fingers and y represents the distance from the backbone. 

A suitable definition of the finite Van Hove sphere for the comb lattice is the square periodic comb of L x L 
sites. In this case the backbone and each finger are rings of L sites. Each vertex can be labelled by (x,y), x, y S 
[0, 1, . . . , L — 1,L], with the conditions (0, 0) = (L, 0) and (x, 0) = (x, L) in order to guarantee periodic boundaries. 

In the thermodynamic limit the spectrum of the infinite comb is evaluated as the limit for L — > oo of the spectrum 
of the hopping model on the finite L x L combs. The eigenvalue equation ( |Io[ ) on the finite comb with hij = —tAij 
with the previous labeling of sites reads: 

L-l 

-t X! [(<^,x'+i +5 x , x i-i)5 yfl 8o,y' + (<W+i + 6y, y >-i)6 x , x >] ip(x',y') = Eip(x,y) (23) 

x' ,y'—0 

By exploiting the translation invariance in the direction of the backbone, a Fourier transform in the variable x 
reduces (|23J) to a 1-dimensional eigenvalue problem. Let us define: 

*P(k,y) = J2e lkx ^(x,y) (24) 

X 

with k — 2nn/L, n = 1 . . .L. The eigenvalues equation ( |23|) becomes: 

- t [-2 COs(k)S y fi6oySkM' + {5y,y'+l + 5y, y >-i)h,k>] ip(k' , y') = Eip(k, y) (25) 

k',y' 



Now (25) is diagonal in the variable k and can be written as: 



[-2COs(k )5y t0 S ,y> + (Sy,y> + 1 + <W"l)] ^(v') = ^(v) ( 26 ) 

y' 

where ip(k,y) = 5(k — ko)ip{y) with fco = 27rn/L, n = . . . L. Equation ( p6|) can be regarded as a 1-dimensional 
problem of a quantum particle interacting with a potential in the origin, V(ko) = —2tcos(ko). The eigenvalues and 
eigenvectors of (^6|) are obtained by imposing the matching condition in y = for the free particle solutions on the 
chain. Indeed, for y ^= 0, ( p6| ) the wave function satisfying ( p6| ) in y ^= 0, L — 1 are: 

i/j(y) = cos(hy + a) for -2t < E = -2tcos(h) < 2t; (27) 
i/,(y) = Ae hy + Be- hy for E = -t(e h + e~ /l ) < 2t; (28) 
V>(y) = A(-l) y e hy + B(-l) y e- hy for E = t(e h + e - h ) > 2t. (29) 

Imposing ipiu) to be a solution in y = and y = N — 1 introduces restrictions on the parameters h, a, A and B, and 
on the eigenvalues E. 

Let us first consider the conditions for (p7|), in y = L — 1 and y = 0: 

-tcos[h(L -2) + a] -icos(a) = -2t cos(h) cos[h(L - 1) + a] (30) 
— t cos[h(L — 1) + a] — tcos(h + a) — 2tcos(fco) cos(a) = —2tcos(h) cos[h(L — 1) + a] (31) 

The system has odd solutions with a = n/2 (i.e. ?p(i) = sin(/?i)), h = 2i:n/L and n — 1, . . . ,L/2 — 1 and even 
solutions, obtained by solving: 



G 



- cos(fco) cot(ft£/2) = sin(ft) (32) 

Equation ([32|) can be solved graphically (see figure 2) obtaining L/2 solutions. In the large L limit the allowed values 
for h are: h rs 7r(2n — 1)/L with n = 1, . . . , L/2. Each values of fco corresponds to L— 1 solutions with energy between 
-2t and 2t 

Let us now consider cases (^,^). For E < 2t the wave function in y = L — 1 and y — satisfies the conditions: 

- ^Ae-' 1 ^- 2 ' + Be h< - L -V) - t(A + B) = -t(e h + e"' 1 ) (Ae"^- 1 ) + Be h ^-^) (33) 
-tiAe-^-V + Be' 1 ^-^) - t(Ae-' 1 + Be' 1 ) - 2 cos(k )(A + B) = -t(e k + e- k )(A + B) (34) 



leading to: 



cos(fco) coth(ft£/2) = smh(h) (35) 



Equation d35| ) can be solved graphically. For cos(fco) > 0, (|35 | ) has a real solution. In the limit L — ► oo this is 
sinh(/i) = cos(fco), corresponding to energy is E = — 2tyl + cos 2 (fc ). If cos(fco) < 0, only the solution with E > 2t 
is present and for cos(fco) = one has the constant solution ip(y) — 1 with energy E = 0. 

The complete spectrum of the comb is now obtained by considering the eigenvalues of (|26| ) for the L values of fco. For 
each fc there are L — 1 eigenvalues of type ( pTj ) with energy 23 = — 2t cosjh) and wave functions ip(x, y) — e lk ° x sin(hx) 
and ip(x, y) = e tkax cos(h\y\ + a), with h and a satisfying conditions ( pfj| , ^ll) . The fraction of states in this spectral 
region is / = L(L — 1)/L 2 . Notice that / tends to 1 in the thermodynamic limit. The normalized density of states in 
this spectral region is: 

dn _ Ld(Lh/2n) _ dE 

with E G [—2t, 2i\. Since the fraction of states with energies E < — 2t or E > 2t vanishes in the thermodynamic limit. 
( |36| ) represents the normalized density of states p(E) of the pure hopping model on the comb graph. 

However, the lowest energy state is obtained from a solution of type (pq). These solutions appear only for cos(fc ) > 
and they satisfy relation (|35|). Since the energy is a decreasing function of h, the lowest energy level is obtained for 
cos(fco) = 1. In this case we have coth(/iL/2) = sinh(/i), and in the thermodynamic limit this corresponds to 
sinh(/i) = 1. From (p8|) we get E — —ty/8. The lowest energy in the normalized density of states is E m = —2t 
and therefore E < E m . Then the eigenvalue equation ( p3| ) presents a hidden spectrum in the thermodynamic limit, 
and the pure hopping boson model on the comb lattice exhibits BEC at low enough temperatures. From equations 
(^8|,[53],[34]) it is also possible to show that on the infinite graph the eigenvector corresponding to the lowest energy 
eigenvalue is ip{x, y) — e~' i ' y '. The wave function for the condensate is localized along the backbone and it decreases 
exponentially along the finger (see figure 3). 

This model does not present an energy gap between Eg and E m , contrary to the typical instance of non interacting 
bosons trapped in a harmonic well. Indeed for each value of fco (cos(fco) > 0) there is a solution of ([35]) with a different 
energy in the interval [Eo,E m ]. In a finite comb of L 2 sites there are L/2 solution of this type and for L — > oo these 
solutions fill densely the interval [E^E m ]. Normalizing the density of states to L, we can obtain the spectral density 
in the thermodynamic limit for E G [EoE m ]: 

dn = d{Lk /27T) = \E\dE 
L L 2vrV8t 2 - E 2 VE 2 - At 2 

An analogous equation holds for the spectral region at high energy (E = [2t,y/St]), where the other hidden states 
appear. The density of states can then be represented as in figure 4, where we used the two different normalizations 
for the hidden spectra and for p(E). 

^,From the knowledge of the complete spectral density of the comb, the thermodynamic quantities relative to BEC 
such as the filling of the ground state as a function of T, the specific heat and the critical temperature as a function 
of / can be analytically calculated . 

V. BOSE-EINSTEIN CONDENSATION ON GRAPHS: FURTHER EXAMPLES 

The spectrum of the pure hopping model can be studied in detail in a wide class of fractals and inhomogeneous 
discrete structures. In particular, the behaviour observed in the comb lattice is typical of bundle structures p0|. This 
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class of graphs are obtained by a "fibering" procedure, i.e. attaching a fiber graph to every point of a "base" graph. 
They are characterized by a hidden low-energy region, and from the general result of section III they exhibit BEC 
at finite temperature. Moreover, as in the comb graph, the wave function of the condensate is localized along the 
base and presents a fast decay along the fibers. An example of this behaviour is found in the brush graph, which is 
illustrated in figure 5. 

An interesting example of a graph where the geometrical setting gives rise to an isolated eigenstate in the low-energy 
region of the pure hopping model is the star graph. The effect, which is analogous to that caused by an impurity 
on a line, here is not produced by the presence of an external potential but simply by the topology of the system. 
Interestingly, the wave function is here completely localized in the center of the star, with an exponential decay along 
the arms. In figure 6 the star graph with the ground state and the spectrum of the pure hopping model are shown. 

Graphs with constant coordination number are typical examples in which the pure hopping model do not present 
hidden low energy regions. This is due to the fact that for this class of graphs the spectrum of the model can be 
obtained from that of the Laplacian matrix by a shift of the zero in the energy. The existence of BEC is therefore 
determined by the spectral dimension d = a of the graph, describing the density of eigenvalues of the Laplacian at 
low energy. This parameter can be exactly calculated for a wide class of discrete structures. On lattices, d coincides 
with the usual Euclidean dimension and one recovers the classical result for BEC on translation invariant structures. 
For exactly decimable graphs (e.g. the Sierpinski gasket Jl^ and the T-fractal, see figure 7) where one always has 
d < 2 there is no Bose-Einstein condensation. 

A fundamental property of the spectral dimension is its independence from the local details of the graph, i.e. 
universality The value of d and the behaviour of the density of eigenvalues in the low energy region is not 
changed under a wide class of transformations, called isospectralities, which can strongly modify the geometry of the 
graph. A simple consequence of this property is that if we consider the pure hopping model on a graph with constant 
coordination number, which differs from a graph of known dimension d by an isospectrality, BEC occurs only if d > 2. 
An example of this behaviour is given by the ladder graph (figure 8) which can be obtained from the linear chain by 
the addition of finite-range links. On this structure then the pure hopping model does not exhibit BEC. 

An important properties of d is that the dimension of the graph obtained as a direct product of two graphs is the 
sum of the dimensions of the original structures. An example is illustrated in figure 9, where we show the direct 
product of a linear chain and a well known fractal, a Sierpinski gasket. In this the coordination number case Zi 
is constant, and applying the previous properties it is easy to show that d = 1 + 21n(3)/ln(4) > 2, and therefore, 
applying the previous results, one immediately infers that BEC occurs on this graph. 



VI. CONCLUDING REMARKS 



We studied in detail BEC of non-interacting bosons on general discrete networks, i. e. on graphs. We evidenced 
that these systems support spatial condensation of bosons into a single state even if d < 2. 

Due to recent advances in engineering technology, Josephson junction networks can be build in a variety of con- 
trollable, non-conventional architectures. The theoretical interest for these systems stems naturally from the ability 
of experimentalists in varying the properties of the networks by acting directly on their geometry. We expect that 
non universal quantities - such as the critical temperature for condensation - will depend crucially on the geometry of 
a graph, thus leading to the possibility of engineering network geometries suitable to observe spatial BEC of Cooper 
pairs. 

Our analysis shows that classical inhomogeneous Josephson junction networks build on graphs may support a spatial 
BEC of Cooper pairs induced only by the geometry of the ambient space on which the junctions lie. We feel our 
prediction of the existence of spatial BEC in Josephson junction networks may be brought soon to direct experimental 
testing. 

Acknowledgements We thank F. Illuminati and A. Trombettoni for enlightening discussions and a careful reading 
of the manuscript. We thank F. Illuminati for suggesting the star graph as an interesting geometrical structure. 



[1] Dalfovo F, Giorgini S, Pitaevskii LP and Stringari S, 1999 Rev. Mod. Phys. 71 463 

[2] Huang K, Statistical Mechanics 1963 Wiley, New York 

[3] Alexander S and Orbach R 1987 J. Phys. (Pans) Lett. 92 108 

[4] Burioni R and Cassi D 1996 Phys. Rev. Lett. 76 1091; 1997 Mod. Phys. Lett. B 11 1095 



8 



[5] For a review see Nakayama T, Yakubo K and Orbach RL 1994 Rev. Mod. Phys 66 381 
[6] Hattori K, Hattori T and Watanabe H 1987 Prog. Theor. Phys. Suppl. 92 108 

[7] Burioni R, Cassi D, Meccoli I, Rasetti M, Regina S , Sodano P and Vezzani A 2000 Europhys. Lett. 52 251 
[8] Burioni R, Cassi D, Vezzani A, |cond-mat 0104365 

[9] Barone A and Paterno G, 1982 Physics and applications of Josephson effect, Wiley, New York 
[10] Simanek E, 1994 Inhomogeneous superconductors: Granular and quantum effects, Oxford University Press 
[11] Fisher MPA, Weichman PB, Grinstein G and Fisher DS 1989 Phys. Rev. B 40 546; Bruder C, Fazio R, Kampf A, van 
Otterlo A and Schon G 1992 Phys. Set. 4 2 159 Scalettar RT , Batrouni GG, Kampf AP and Zimanyi GT 1995 Phys. Rev. 
B 50 8467; Fazio R and van der Zant H, |cond-mat/0001152 
[12] van Oudenaarden A and Mooij JE 1996 Phys. Rev. Lett. 76 4947 

[13] Jaksch D, Bruder C, Cirac JI, Gardiner CW and Zoller P 1998 Phys. Rev. Lett. 81 3108; Trombettoni A, Ph. D. Thesis, 
SISSA (2001) 

[14] Trombettoni A and Smerzi A 2001 Phys. Rev. Lett. 86 2353 

[15] Anderson PB and Kasevich MA 1998 Science 282 1686; Cataliotti FS, Burger S, Fort C, Maddaloni P, Minardi F, 

Trombettoni A, Smerzi A and Inguscio M 2001 Science 293 843 
[16] Harari F 1969 Graph Theory (Addison- Wesley Reading); Biggs NL 1974 Algebraic Graph Theory (Cambridge University 

Press); Mohar B and Woess W 1989 Bull. Lond. Soc. 21 209 
[17] Burioni R, Cassi D and Vezzani A 2000 Eur. Phys. J.B. 15 665 
[18] Burioni R, Cassi D and Vezzani A 1999 Phys. Rev. E. 60 1500 
[19] Rammal R 1984 J. Physique. 45 191 
[20] Cassi D and Regina S 1997 Mod. Phys. Lett. B 11 997 



FIG. 1. The comb lattice. 




FIG. 2. The numerical solution of equation (|32j) for ko = and L = 16. We represent the function cot(L/i/2) with a 
continuous line and s'm(h) with a dashed line. The intersections of the two curves are the L/2 solutions of (133). 







FIG. 3. The wave function of the low energy eigenstate in the comb lattice. 



FIG. 4. The spectrum of the pure hopping model on the comb graph. p(E) (normalized to L 2 ) is plotted with a continuous 
line, while for the hidden spectra (normalized to L) we used the dashed line. 




FIG. 5. The brush graph and the spectrum of the pure hopping model. As in the comb graph, different normalizations are 
used for the hidden spectra (dotted line) and p(E) (continuous line). 
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FIG. 8. The ladder graph and the spectrum of the pure hopping model on this structure. Here we do not have any hidden 
spectrum and p(E) — > oo when E — > E m = Eq (d = 1). 
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FIG. 9. The graph obtained as a product of the Sierpinski gasket and a linear chain. Since the coordination number is 
constant, there is not hidden spectra in the low energy region. However d — 1 + 2 ln(3) / ln(4) > 2 (sum of the dimension of the 
original graphs) and therefore BEC occurs on this structure. 
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